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ABSTRACT 

Although OGLE-TR-56b was the second transiting exoplanet discovered, only one light curve, 
observed in 2006, has been published besides the discovery data. We present twenty-one light curves 
of nineteen different transits observed between July 2003 and July 2009 with the Magellan Telescopes 
and Gemini South. The combined analysis of the new light curves confirms a slightly inflated planetary 
radius relative to model predictions, with R p = 1.378 ± 0.090 Rj. However, the values found for the 
transit duration, semimajor axis, and inclination values differ significantly from the previous result, 
likely due to systematic errors. The new semimajor axis and inclination, a — 0.01942 ± 0.00015 AU 
and i — 73.72 ± 0.18°, are smaller than previously reported, while the total duration, T u — 7931 ± 38 
s, is 18 minutes longer. The transit midtimes have errors from 23 s to several minutes, and no evidence 
is seen for transit midtime or duration variations. Similarly, no change is seen in the orbital period, 
implying a nominal stellar tidal decay factor of = 10 7 , with a three-sigma lower limit of 10 5 7 . 
Subject headings: stars: planetary systems - OGLE-TR-56 



1. INTRODUCTION 



OGLE-TR-56b was the se cond transiting exoplanet R and V filters 1 Pont et a M 2007 D 



discovered, after HD209458b (Charbonneau et al. 



single light curv e in 2006 observe d alternately in Bessel 

The precision of the 



2000 



Henry et al.||2000 1 and the first found by a photometric 



wide-field survey (Udalski et al. 2002). The planetary 



nature of this object was confirmed by |Konacki et al.| 



(2003), who derived a mass of 1.3Mj from radial veloc- 
ity measurements of the host star. An initial flurry of 
follow-up observations i nclud ed ra dial velocity measure 



planetary parameters for this system has thus far been 
limited by the p hotometric quality of the available tran- 
sit light curves ( Southworth |2008 ) and not by the preci- 
sion of the stellar parameters, as is currently the case for 
many other planets. 

OGLE-TR-56b was selected in 2006 as part of our cam- 
paign of high-quality photometric observations to search 



merits by Torres et al. (2004) and Bouchy et al. (2005); , , 

the deter mination of th e~Tioit star' s fundamen tal rJSriS- |et al.||2010a|b| [2011] . Its short orbital period allows tor 



for transit timing v ariations (see Adams 2010 Adams 



2006); 



eters and chemical composition by Santos et al 
and more re cently the detection of OGLE-T R-56b's at 
mosphere by Sing and Lopez-Morales (2009). However, 
due to the faintness of the target {V = 16.56, / = 15.3) 
and the recent explosion in the number of other, brighter 
transiting systems, OGLE-TR-56b remains relatively un- 
studied, despite having one of the shortest known periods 
(1.2 days). The radius of the planet is currently con- 
strained by only two light c urves: the composite OGLE 
discovery light curve in / ( |Udalski et al.||2002[ ), and a 
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many observational opportunities, and it remains one of 
the better candidates to search for changes in the orbital 
period due to orbital decay. Predictions vary widely on 
the time scale over which orbital decay would be ob- 
servable (|Sasselov| 2003: Piitzol d^et al.p 
Patzold|2007||Levrard et a l.|2009p , with 



2004 Carone and 



estimates of the 

remaining lifetime o f this system ranging from a few bil- 
lion years (Sasseloy 2003) t o a much shorter 7 Myr (as 
ird 



suming = 10° |Levrard et al. 
assumptions about tidal equili 



2009 1 using very different 
Dnum modes. Therefore, 
the estimates for the decrease in the orbital period due 
to orbital decay range from 0.1 to 10 ms yr _1 . This value 
is only a few times smaller than the current error on the 



orbital period (86 ms fromlPont et al. 



2007) 



Here we present twenty-one new light curves of nine- 
teen transits of OGLE-TR-56b, observed from 2003 to 
2009, which we use to improve the planetary system pa- 
rameters and also to search for potential timing anoma- 
lies in this system, including transit timing variations 
(TTVs) caused by interactions other bodies in the sys- 
tem, and transit duration variations (TDVs) caused by 
orbital precession, among other effects. In § [2] we de- 
scribe the collection and analysis of the data. In § [3] we 
describe the transit model fitting; we note the difference 
in light c urve shape compa red t o th e light curve pub- 



lished by Pont et al. (2007) in § 3.3 In §[4] we present 



the revised transit ephemeris and discuss how the timing 
constrains the presence of other objects in the system. 
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2. OBSERVATIONS AND DATA ANALYSIS 

Nineteen transits were observed between 2003 July and 
2009 July, all but one with the Magellan telescopes. All 
transits are referred to by the UTC date of observations; 
when multiple light curves exist for a given UTC date, 
e.g. the two light curves on 2005 April 19 UT in the 
Johnson B and I band, they are denoted respectively as 
20050419B and 200504191. Seven full or partial transits 
were observ ed by the Transit Lig ht Curve (TLC) collabo- 
ration (e.g., Holman et al.|2006 ) during 2003-2006, using 



MaglC-SITe. Thirteen transits were observed between 
2006-2009 on Magellan using three instruments: POETS 
(3 transits), IMACS (1 transit), and MagIC-e2v (9 tran- 
sits). One additional transit was observed in 2009 using 
CMOS on Gemini South. Details on instrument setup 
and observational conditions are summarized in Table [TJ 
The Magellan instruments all had relatively small 
fields-of-view with high spatial sampling. The earliest 
transits used the Magellan Instant Camera (MagIC) with 
the original SITe CCD on the Clay telescope. MaglC- 
SiTe had a field of view of 140" x 140" and a plate scale 
of 0.069" per pixel; the CCD gain was about 2 e-/ADU 
and the read noise 6 e-. Four filters were used: Harris 
B, Mould /, and Sloan r' and z' (the filter transmission 
curves are shown in the legend of Figure [TJ). In June and 
July 2006, POETS (Portable Occultation, Eclipse, and 
Transit System) was installed as a visiting PI instrument 
on the Clay telescope. POETS has a field of view of 23" 
x 23" and a plate scale of 0."046 per pixel. The camera 
is a frame-tra nsfer CCD which ca n be GPS triggered and 
is described in 



Souza et al. 



( 2006| ). POETS was operated 
full-frame (no binning) in conventional mode with a gain 
of 3.4 e-/ADU and read noise of 6 e-. For the first two 
transits observed, a Schuler Astrodon Johnson-Cousins 
Is filter was used. For the third transit, observed on 
20 Jul y 2006, coincident ally the same night as the tran- 
sit by Pont et al. (2007), no filter was used as a test to 
achieve greater throughput; the POE TS CCD is sensitiv e 
to wavelengths between 400-900 nm ( |Gulbis et al.|2008[ ). 
We refer to this a s 20060720N to distinguish it from the 
Pont et al. (20071 observations. 

One transit was observed in 2007 with the Inamori 
Magellan Areal Camera and S pectrograph (IMAC S) on 
the Magellan Baade telescope ( |Dressler et al.|2~006 ). We 
used the //4 imaging mode with a subraster on one of the 
eight instrument CCDs (chip 2), reading out 1200x1200 
pixels centered around x=1080 and y=3015, in order 
to decrease the readout time to roughly 35 seconds per 
frame in FAST mode. The area imaged was thus 133" 
x 133", with a plate scale of 0.11" per pixel. This chip 
has a gain of 0.9 e-/ADU and read noise of 4.6 e-. A 
wideband filter, WB6300-9500, was used. 

Between January 2008 and December 2009, the Magel- 
lan Instant Camera (MagIC) was mounted on the Magel- 
lan Baade telescope with the addition of a new e2v CCD. 
The MagIC-e2v detector, which shares a dewar with the 
older SITe CCD, is identical to the red CCD on HIPO 
(one of the first generation instruments to be flown on 
SOFIA), a fast read-out dir ect imaging camera that uses 



the LOIS control software (|Dunham et al. 2004 Taylor 



et al.|2004H"Osip et al.|2008l 

sits were observed in the Sl( 



of 0."037 per pixel unbinned. The camera can be oper- 
ated in two different modes: single exposure mode, with 
a readout time of about 5 seconds per exposure, which 
was used for the first three transits in 2008; and frame 
transfer mode, with a readout time of only 3 millisec- 
onds between frames in an image cube, which was used 
for the following six transits. The gain and read noise of 
the transits observed in 2008 were 2.4 e-/ADU and 5.5 
e- per pixel, respectively, and two amplifiers were used 
during readout; after July 2008, the CCD was reconfig- 
ured to have a gain of 0.54 e-/ADU and 5 e- read noise 
per pixel and a single readout amplifier. 

The transits observed with the frame transfer camera 
POETS were observed at a very rapid cadence (2-6 s per 
exposure), in order to get good time sampling during 
transit ingress and egress. However, the noise per frame 
was such that we chose to bin every 10 frames (20 s) for 
20060622, 8 frames (32 s) for 20060714, and 10 frames (60 
s) for 20060720N. In subsequent transits, the exposure 
times were adjusted to maintain roughly 10 6 integrated 
photons for both the target and multiple nearby com- 
parison stars. Exposure times for unbinned (1 x 1) data 
ranged from 17-60 sec, and for binned data from 10-60 
sec. Details of the observing settings are noted in Ta- 
ble [l] Gaps in the data for transits 20060622, 20060714, 
and 20080727 were caused by separate instrument com- 
puter glitches, while the gap in transit 20090510 was due 
to observing a gamma ray burst for another project while 
OGLE-TR-56b passed through zenith. 

The observation strategy for the transits observed by 
the Transit Light Curve collaboration between 2003-2006 
is simila r to that of OGLE-TR-lllb. Details are pro- 
vided in 

was observed using GMOS 



Winn et al.| (2007 
The transit on 2009TJ42" 



All of the MaglC-e2v tran- 
oan i' filter. The MagIC-e2v 
camera has a field of view of 38" x 38" and a plate scale 



(Gemini Multi-Object Spectrograph) at Gemini South 
Telescope. The observations were done in service mode 
using a Sloan i filter (LG0327), 2 x 2 binning, and an 
integration time between 12 and 30 seconds. We used 
the Region of Interest (ROI) of the camera in order to 
reduce the readout time to only 24 seconds in fast readout 
mode and low gain. 309 frames of 1024x2304 pixels were 
obtained with a scale of 0."146 per pixel. The chip has 
a gain of 5 e-/ADU and read noise of 7.8 e-. 

Accurate timing is of the utmost importance for this 
project, so special care was taken to ensure that the 
correct times were recorded in the image headers. For 
the transits observed with POETS in 2006, each image 
frame was triggered by a GPS, so the UTC start times 
are accurate to the microsecond level. For the transits 
observed with GMOS, MaglC-SITe, IMACS and MaglC- 
e2v (in 2008), the UTC start times for each image were 
recorded from network time servers, which in most cases 
were verified by eye to be synchronized with the observa- 
tory's GPS clocks at the beginning of each night. For the 
MagIC-e2v observations in 2009, the times came from a 
small embedded control computer (a PC 104), which re- 
ceived unlabeled GPS pulses every second. As with the 
network time servers, the PC104 was synchronized with 
the observatory's GPS before each transit observation. 
In all cases the UTC time signals written to the image 
headers agree within one second with the GPS time, as 
verified by examining the system control logs. 

2.1. Data analysis 
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All data were calibrated using IRAF 10 . The zero level 
for POETS and IMACS data frames were calculated from 
bias frames taken before or after the transits, while the 
MagIC data frames were corrected using overscan regions 
on each image. The images were flat-fielded using either 
dome or twilight flats in the appropriate filter, as avail- 
able. The CMOS data were reduced using the Gemini 
pipeline. 

Although OGLE-TR-56 is in a very crowded field (Fig- 
ure [5]) , the generally excellent seeing and good spatial 
resolution of the detectors allowed for high-quality light 
curves using aperture photometry. Most of the photom- 
etry was done using the IRAF routine phot, part of the 
apphot package. A wide range of apertures and differ- 
ent comparison stars were examined, depending on the 
nightly conditions (e.g., seeing) and specific field of view, 
in the sa me way as described for othe r planets we have 
analyzed ( |Adams et al.||2010a|b[ |2011| ). 

For a few transits, an alternative aperture photometry 
method was used to achieve greater precision. In this 
method, boxes were drawn around the target and com- 
parison stars, and the sky was determined using a single 
30-pixel box drawn around an empty region of sky (not 
easy to find for targets in crowded fields, such as OGLE- 
TR-56b). This method, as implemented in Mathematica, 
is both slower and not as robust as IRAF's phot, partic- 
ularly for data sets in which the stars shifted by more 
than a few pixels; but for two transits (20060622 and 
20080514) it produced lower out-of-target scatter and 
cleaner light curves. The sky region for 20060622 was 
centered 2.1"west and 0.1"south of the transiting planet 
host star; for 20080514 the region was located 2.0"east 
and 0.3"south. 

We additionally explored using image subtraction 
( |Alard and Lupton||1998{ |Alard|2000| ) to create the light 
curve for one transit, 20090504. This test produced a 
light curve with slightly better scatter compared to the 
aperture light curve (1.0 vs 1.1 mmag in 2 min), but 
the depth of the image-subtracted light curve is much 
shallower (k = 0.092 ± 0.002) than the aperture light 
curve for either the individual fit for that light curve 
(fc = 0.102 ±0.004) or the joint fit (k = 0.1039 ±0.0004). 
This problem is a known potential pitfall fo r image sub- 
traction photometry (see e .g. discussion in Gillon et al 



2007 Adams et al. 2010a). We note, however, that im 



age subtraction does not always give erroneous depths; 
the TLC light curves w ere created with image subtrac- 
tion (Winn et al. 2007), and the depths and other pa- 
rameters, particularly of the full transits, agree well with 
light curves made with aperture photometry. 

After obtaining light curves for each transit, the out-of- 
transit flux of each curve was examined for linear trends 
against several variables: airmass, seeing, telescope az- 
imuth, (x, y) pixel location, and time since beginning of 
transit. These parameters were chosen because they are 
either directly correlated with photometric trends (e.g. 
seeing, airmass), or are proxies for other effects that may 
be harder to measure (e.g. the telescope azimuth, which 
is a major component in the de-rotator rates, which were 

10 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



not recorded in the headers of all transit files). The 
data presented no significant systematics, except in three 
transits, for which we fit and remove linear trends for 
slight slopes in the flux over time (transits 20070830 and 
20090504), and a slope with seeing (transit 20080727). 
Since these slopes were removed before the light curve 
parameter fitting step (see § k3|) , for simplicity, it is possi- 
ble that covariances between those trends and the transit 
parameters might have a slight effect on the individual 
transit parameter. However, the joint fit values are un- 
likely to be affected since most of the transits fit had no 
slopes removed. 

In the case of 20060720N, the light curve observed with 
no filter, we found that differential color extinction domi- 
nates the light curve systematics. Differential extinction, 
strongly correlated with airmass, is particularly trouble- 
some, since it affects individual comparison stars differ- 
ently We therefore took a slightly different detrending 
approach for this dataset. The low airmass values (1.0 
to 1.1) allow a linear approximation of the airmass term 
when detranding. In addition, we fit a separate linear 
slope to the ratio of the target to each comparison star, 
before attempting to combine them to produce the final 
transit light curves. Although we initially examined 13 
different comparison stars over a wide range of photomet- 
ric apertures, only two comparison stars could be used in 
the final light curve. We note that despite the relatively 
low scatter of the resulting light curve, there are dangers 
in this approach, notably in how to accurately detrend 
when there is little data on one side of the out-of-transit 
baseline to constrain the slope removed. However, given 
that the resulting light curve agrees well in depth and 
shape with the other, filtered transits, this method seems 
to have worked in in this case. Because correcting for dif- 
ferential extinction in bluer wavelengths introduces many 
systematic complications, we do not recommend observ- 
ing planetary transits without filters as a general rule; if 
the goal is to achieve greater throughput, a redder wide- 
band filter would be a better option. 

The photometric parameters used to generate each 
light curve are shown in Table [2j Each resulting light 
curve is plotted in Figure [TJ together with the residu- 
als from the best model plot, with each transit binned 
to two minutes for comparison. The transit fluxes and 
times are available as an online table; an excerpt is shown 
in Table H 

2.1.1. Literature light curves 

The only published high-quality light curve for OGLE- 
TR-56b was observed on 2006 July 20 UT with the VLT 
in both R and V , with alternatin g sequences in ea ch filter 
of 7-8 exposures of 15-40 s each ( Pont et al.|2007[ ). Upon 
refitting that light curve (provided by F. Pont 2007, per- 
sonal communication), we found that the shape of the 
transit does not agree well with the one obtained from 
the new light curves, most notably in an 1 8-minute dis- 
crepancy in the transit duration (see 



3.3) 



In an attempt to resolve this problem, we downloaded 
the now publicly available on the VLT archive database 
raw image frames and re-analyzed the data for each filter 
(see Figures 1 and 3). Our new light curves, based on 
aperture photometry, are less precise than the original 
image subtraction light curve: 1.0 mmag in i?-band and 
1.2 mmag in U-band, compared to 0.7 mmag for the orig- 
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inal combined V/R curve. However, the aperture light 
curves are consistent in transit depth and duration with 
the 19 other light curves presented in this paper. The 
revised curves, referred to as 20060720R and 20060720V 
henceforth, are used in all subsequent analyses, except 
as noted for comparison with the original photometry, 
which we call 20060720P. We also note that the times 



Carlo (MCMC) method, us ing Gibbs sampling and 



used in the Pont et al. (2007) light curve appear to be 
the start time, rather than the mid-exposure time of each 
frame, and thus are earlier than the corresponding points 
in our analysis by 8-20 s; however, this slight time offset 
is less than the midtime error on the transit, and cannot 
account for the duration difference. 

We did not refit the OGLE survey light curve, which 
is a composite of 13 fu ll and partial transi ts spanning 



several hundred nights (Udalski et al. 20021). However, 
we do use the most up- to-date pu blished nrldtime of that 
composite light curve (Torres et al. 20041 in our timing 
analysis, after adding 66.184 s to correct for the UTC-TT 
offset (§[§. 

3. TRANSIT FITTING RESULTS 
3.1. Model 



Our trans it model fits use the algorithm of Mand eT~and| 
Agol| ( |2002'| ) as implemented by |Carter and Winn| ( |2009p , 
assuming white noise. We assumed that OGLE-TR-56b 
has zero obliquity, oblateness and orbital eccentricity. 
The stellar mass and radius values used to convert the 
model output into p hysical parameters were taken from 
Torres et al. ( 2008 ) and are listed in Table |4l We also 



fixed the orbital period to P = 1.211909 daysl fPont et al. 
2007 1 , since that parameter is only used to convert a/R* 



into an orbital speed, and therefore has no major effect 
in the other model parameters. To model the stellar 
limb darkening, we assumed a quadratic law with initial 
values for the u\ and ui param eters set to the AT LAS 



2000 2004 1 11 In 



values for the appropriate filter ( Claret 
all transits we found it necessary to fix the quadratic 
term U2- In addition, we only fit for the linear term u\ 
in the transits observed in the i' and /-band filters, while 
leaving it fixed for the transits observed in the B, r' , z' , 
W B-han&s and with no filter because the precision of 
these light curves is insufficient to constrain the value 
of Mi. The values for ui and u? are calculated using 
thcjktld program by Southworth (2008) 12 , assuming the 

stellar parameter values T = 6119 K, logg — 4.2 cm/s 2 , 
[M/H] = 0, and V m i cro = 2 km/s. We fixed the coeffi- 
cients for the wideband transit 20070830 to be the same 
as the Sloan i! filter, since the actual filter used, WB6300- 
9500, is centered at the same wavelength, though about 
twice as wide. The limb darkening parameters for the 
transit with no filter were derived by setting all other 
parameters fixed to a joint-fit value excluding this light 
curve, then fitting for ui tTione and U2,none alone. 

3.2. Light curve fits 

Each light curve was fit both independently and jointly 
with all other light curves using a Markov chain Monte 

11 Initial values: ui none = 0.19, «2 none = 0.25, u\ b = 0.5458, 
«2, s = 0.2504, u iy ' = 0.4089, u 2 y = 0.2757, u iy = 0.2289, 
«2, z ' = 0.3112, u\j = 0.1958, «2j = 0.3561, u-^y = ui^wb = 
0.2146, u 2 i> = U2 wb = 0.3569, «i v+R = 0.3157, J12 v+R = 
0.35265. 

12 http://www.astro.keele.ac.uk/ jkt/codcs/jktld.html 



Metropolis-Hasti ngs stepping ( Tegmark et al.||2004 Hol- 
man et al. |2006[). This me thod is described m greater 



dot a: 



11 m e.g. 



Uarter et al. (20111. Three independent 



chains of a million links each were combined, discarding 
the first 50,000 links of each chain, to avoid any potential 
solution biases due to the initial values of the input pa- 
rameters. For each fitted parameter - the radius ratio, k, 
inclination, i, semimajor axis ratio, a/ R„ (global param- 
eters for all transits), limb darkening coefficients u\ and 
U2 (global for each filter), and out-of-transit flux, Foot 
and transit midtime, Tq (specific to each transit) - we 
find the posterior distribution, all of which are roughly 
Gaussian. The median value and 68.3% confidence in- 
tervals of each parameter distribution are reported in 
Table H 

To account for excess correlated noise in the light 
curves, we calculated the time-averaged residuals. This 
is done by binning the residuals for each light curve into 
bins from 10 to 30 minutes, the typical time scales of 
correlated noise (Pont et al.||2006[), and then calculating 



how much greater the actual noise is compared to the 
ideal noise assuming photon statistics. The excess noise 
varied from 1-2.8 times the predicted noise level, and has 
been included in the error bars on all transit midtime and 
individual parameters. 

As a check that there were no additional slopes in the 
data that could affect the photometry, we ran a test 
MCMC fit to 12 light curves that included a slope with 
respect to airmass as an additional free parameter. None 
of the 12 full light curves fit in this way had a significant 
shift in transit parameters, so we did not included air- 
mass slopes in the rest of our analysis. 

Table [4] summarizes the joint model fit values for the 
planet-star radius ratio, k, inclination, z, semimajor axis 
ratio, a/i?*, total transit duration, T14, full transit dura- 
tion, T23, and ingress/egress duration, T 12 = T34. (The 
total transit duration is the time from first to last con- 
tact, while the full transit duration is the time when the 
disk of the planet is entirely overlapping the disc of the 
star. The ingress and egress duration are the times it 
takes for t he transition and are equal if the orbit is cir- 
cular. See Winn (2010) for a mathematical description 
of the transit nomenclature.) 

We also fit each light curve independently. Table [5] 
reports each individual parameter and error, along with 
the time-averaged residual factors by which those errors 
have been increased. To explore variations with time, the 
values are plotted in Figure [3J with the joint-model fit 
shown as a horizontal line with l-cr errors. A few partial 
or noisy light curves (20050419B, 200504191, 20060714) 
have been omitted from the individual results; similarly, 
the radius ratio for 20030730z was fixed to the joint-fit 
value. In general, the independent fit resultss are con- 
sistent with the joint fit, with the notable exception of 
the original photometry for 20060720P, shown as a red 
triangle. No significant variations are seen in any of the 
studied parameters over the 6-year time span of the ob- 
servations. 

3.3. Comparison to previous system parameters 

With twenty-one new light curves, the combined fit 
provides a more precise determination of the system pa- 
rameters, as summarized in Table |4j The errors derived 
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from the photometry for the radius ratio, inclination, 
and semi-major axis are now significantly smaller than 
the error in the measured stellar radius; the component 
of the error in the planetary radius ratio due to stellar 
noise is fifteen times greater than the contribution from 
the photometry. Any improvement in the stellar param- 
eters will significantly improve the precision with which 
the radius of OGLE-TR-56b is now known. 

The radi us ratio rem a ins co nsistent with the value re- 
ported by |Pont et al ] ( |2007[ ), based on a single high- 
precision light curve. Moreover, of the nineteen light 
curves with independently fit radius ratio values in Ta- 
ble [5j none lies more than 1.3 a from the joint-fit value. 
Most of the new light curves have radius ratio precisions 
ranging from 1-4%. The lack of change in the transit 
depth over time indicates that either the star is less ac- 
tive than the few percent level, or that the pattern of 
star spots changes very slowly. No evidence was seen of 
star spot crossings. 

Although our radius ratio agrees with |Pont et al.| 
(2007), the transit duration, T14 is significantly longer, 
by 18 minutes, with corresponding differences also seen 
in the highly-correlated inclination and semi-major axis. 
The duration difference can be seen by eye, as is illus- 
trated in Figure [3] with three light curves of the same 
transit on 2006 July 20. The top curve in Figure [4] cor- 
responds to the original photometry reported in |Pont| 



et al.| ( |2007| ) for the obser vations with the V LT. The sec 
ond curve snows the same Pont et al. (2007) data frames, 
but with new aperture photometry m the same manner 
as the other transits reported in this paper. The third 
curve shows the same transit epoch observed indepen- 
dently on Magellan. The solid (black) and dashed (gray) 
lines in each panel show, respectively, the best fit to all 
new light curves, and the best fit to the |Pont etaT] ( |2007[ ) 
light curve alone. A fourth curve for a different transit 
(20080514) is shown to illustrate the good agreement be- 
tween the aperture photometry light curves and the rest 



of the data, though not the photometry of |Pont et al. 
((2007) 



The duration of all transits, except for 20060720P, are 
best fit by the longer-lasting transit model (black solid 
line) . We th us conclude that the photometry from |Pont 



et al. (2007) most likely suffers from some unidentified 



noise source that shortens the apparent duration of the 
transit; this problem was likely exacerbated by the lack 
of data immediately before transit. This finding should 
be taken as a caution against using single-transit light 
curves to argue for transit duration variations due, for 
example, to planetary orbital precession. 

The new value obtained for the semi-major axis ratio 
of the system, a/Ft* = 3.065 ± 0.022 ve rsus the value of 
3.75 ±0.15 derived by |Pont et aT] ( |2007[ ) , means that the 
planet is somewhat closer to the star than previously es- 
timated. The implications for the planetary temperature 
are discussed in 



3.4 



The new value for the inclination, 
1 = 73.72 ±0.18, also means that the planet is closer to a 
grazing orbit than previously thought (b — 0.859±0.003). 

We also note that the revised photometry should 
prompt a self-consistent reanalysis of the stellar param- 
eters, which are now mildly discrepant. OGLE-TR-56 
is faint (V = 16.56), and the radial velocity m easure 



ments have uncertainties from 20 — 100 m s (Torres 



et al. 2004 Bouchy et al. 2005). Therefore the am- 



plitude of the radial velocity signal, and subsequently 

Torres 



the mass of the star, are not well constrained. 



et al. 



(20081 derived a value of M* = 1.228 ± 0.00 



78Mr 



based on the stellar densit 



derived from the photom- 
etry of Pont et al. (20071). Combined with the or- 
bital period, f = 1.2119819 days, this implies a semi- 



major axis of 3.76 



+0.35 

-0.31' 



inconsistent with our new result 



(a/Ft* = 3.065 ± 0.022). The new photometry suggests a 
lower stellar density is required, with shifts of order —la 
in mass and ±3 a in stellar radius needed to achieve con- 
sistent values for a/R*. A re-analysis of the the radial 
velocity data, perhaps including new measurements, are 
recommended to accurately determine the stellar param- 
eters. 

3.4. Implications for the observed planetary occultation 



Sing and Lopez-Morales ( 2009 ) obtained a z'-band oc- 
cupation depth for OOLJ^TR^56b of 0.0363 ± 0.0091% 
and a brightness temperature of T z > = 27181^7 K, us ~ 
ing the va l ues of the system parameters derived by [Tor- 



res et ah] (2008) and an effective temperature for the 
star of 6119 ± 62K. Since the system parameters have 
changed, we refit the oc cultation data in the same man - 
ner described in §3.1 of Sing and Lopez-Morales (2009), 
replacing the orbital period, a/R*, the planet-to-star ra- 
dius ratio, and the orbital inclination by the new values 
in Table |4j The new fit gives an occultation depth of 
0.0374% and a central phase of = 0.4986 ± 0.0007. 
The central phase is within the errors of the previously 
derived value and still consistent with a circular orbit. 
The occultation is 0.0011% deeper than the value derived 
before, but both values still agree within the measured 
depth uncertainty. The revised brightness temperature 
of the planet is T z > = 2708^^0 K, also consistent with the 
value derived before. Therefor e, the atmospheric prop- 
erties of the planet derived in |Sing and L6pez-Morales| 
(20091 remain valid, unless significantly different values 
for the stellar radius and temperature are found. 

4. TIMING 

With twenty different transit epochs of OGLE-TR-56b 
measured over eight years, the timing of the system can 
be quite well constrained. In Table |6j we present all 
available transit midtimes, in cluding the OGLE survey 
time from |Torres et al. (2004) as well as the new transits 
presented in this work and the re-analyzed data from 



Pont et al.| ( |2007 |. 

We calculate a new transit ephemeris using the tran- 
sit times derived from the joint fit and also including 
the OGLE survey time (transit 20010615). Two tran- 
sits observed in alternating filters (20030730 in r' and z' 
and 20060720 in R and V) had divergent transit mid- 
times when split by filter, so we fit a single combined 
light curve for each transit individually (see also the com- 
bined fit parameter results in Table [5]). In both cases, the 
midtime from the combined light curve lies between the 
times derived from individual filter curves and is closer 
to the expected time of transit; this indicates the poten- 
tial for timing errors if only a partial or poorly-sampled 
light curve is fit. (The individual filter curves for a third 
transit, 20050419, observed in two filters, B and /, was 
left alone, because both filters cover only half a transit 
and yield large overlapping errors; although included in 



G 



the weighted fit, this transit has minimal impact on the 
resulting ephemeris.) 
We find the new transit ephemeris to be: 

T c = 2453936.60070(44) [BJD] + 1.21191096(65) AT. (1) 

where Tc is the predicted centr al time of a transi t (note 
that this is the same epoch as Pont et al. 2007), N is 
the number of periods since the reference midtime, and 
the values in parentheses arc errors on the last digits. 
The residuals from this ephemeris are shown in Figure [5j 
with the lower panel zoomed in on the crowded region 
around the zero epoch. We note that the transit mid- 
time is robustly measured for the zero epoch, with consis- 
tent timing results whether we use the original photom- 
etry (20060720P), the new photometry (20060720V+R) 
for the VLT light curve, or the independent time from 
Magellan (20060720N). The best linear fit has a reduced 
X 2 = 2.6, so we have scaled the errors on the parame- 
ters in Equation [l] by \/2J5 = 1.6. The relatively poor 
reduced x 2 suggests that the errors on some transits are 
still underestimated; removing the most discrepant 3-4 
transits from the fit results in a reduced x 2 °f j us t over 
1. If these errors really are correct, then it is possible 
that the assumption of a constant ephemeris is incor- 
rect. The most deviant points are 20060720N at transit 
number (—3.6 a), and 20090612 at transit number 873 
(—3.4 a), but there is no reason to conclude at this time 
that intrinsic timing variations are responsible. 

4.1. Limits placed on timing variations 

The error on the revised period estimate is 56 ms, with 
no evidence of decrease seen over the eight year time span 
from 2001-2009. If we fit for a linear decrease in the or- 
bital period over time, the rate in change of the period 
is P = —2.9 ± 17 ms yr^ 1 , consistent with a constant 
period. Based on this result, we can derive a conserva- 
tive lower limit estimate for the value of the stellar tidal 
decay factor, Q*. Assuming a three-sigma upper limit 
on the period change (5.4 ms yr" 1 ) and Equation 5 from 
|Levrard et al.|p009[ ), we find that the tidal decay factor 
for OGLE-TR-56 is no less than Q„ = 10 5 7 . Further- 
more, the nominal value obtained for the period change 
implies aQ„ = 10 . 

No evidence is seen of transit duration variations. The 
best fit to the individual transit durations listed in Ta- 
ble [5j using the combined light curves 20030730r+z and 
20060720R+V in place of the curves separated by fil- 
ter, and omitting the much-shorter duration light curve 
20060720P, finds that the total transit duration has 
changed by T\± = —0.3 ± 32 s yr _1 , consistent with no 
change. Since the planet has a large impact parameter 
and nearly-grazing transits, it is possible that in the fu- 
ture orbital precession might lead to observable changes 
in the planetary duration. However, no such changes 
have been seen to date. 

We performed numerical integrations to place limits 
on the perturber mass for a subset of the data, using 
the same method for other planets described in |Adam"s| 
et al.| ( |2010a|bl|2011| . Using 11 transits from 2006-2009 



companion mass limits were placed down to 12Mj in the 
exterior 2:1 mean motion resonance. However, we were 
unable to successfully run the analysis on the full set of 
transit times, most likely owing either to intrinsic insta- 



bility in the system (meaning very few potential compan- 
ions would be stable) or an error in our assumptions (that 
the timing residuals are intrinsically flat). Given the high 
reduced x 2 value of the constant-period fit in Equation!!] 
it is most likely that either a few transit errors remain 
underestimated, or that there are timing variations that 
will require more data to conclusively identify and char- 
acterize. 

5. CONCLUSIONS 

In this work we have presented 21 new light curves 
of 19 transits of OGLE-TR-56b, vastly increasing the 
supply of high quality data on the planet. Our fitted 
radius value of 1.378 ±0.090 Rj is almost identical to the 
previously published value, and we note that the error is 
almost entirely supplied by error on the stellar radius; the 
component of uncertainty from the photometry alone is 
fifteen times smaller. 

The values presented for the transit duration, inclina- 
tion, and semi-major axis in this work are significantly 
different from those reported earlier, the most likely ex- 
planation being an error in the previous photometry. The 
new value for a/R* = 3.065 ± 0.022 places the planet 
slightly closer to its star than previously thought, while 
the inclination, i = 73.72 ± 0.18 degrees, means that 
planet is closer to a grazing orbit (b — 0.859 ± 0.003). 
Furthermore, the photometrically-derived values for the 
semi-major axis ratio are somewhat inconsistent with the 
current values for the stellar radius and mass; the stel- 
lar radius would need to be larger and the stellar mass 
smaller in order for the values to agree. New measure- 
ments and analysis of the stellar parameters are needed 
to resolve this inconsistency. In particular, if the stellar 
radius really is larger, the planetary radius value would 
also increase, implying that OGLE-TR-56b has an even 
more inflated radius than currently thought, with impli- 
cations for its composition and energy budget. 

The new orbital period is P = 1.21191096±0.00000065 
days. Although there are a few points with 3-sigma resid- 
uals, there is no reason to conclude at this time that these 
are transit timing variations, the most likely explana- 
tion being underestimated errors. Excluding the likely 
erroneous low value for the previously published dura- 
tion, the transit duration has remained constant during 
a nine year span (2001-2009), with the best fit change 
T14 = —0.3 ± 32 s yr _1 . The orbital period has likewise 
been constant, with P = — 2.9±17 ms yr _1 , consistent 
with no change. Taking the three-sigma upper limit on 
the period change (5.4 ms yr _1 ) means we can place a 
lower limit on the stellar tidal decay factor of Q* = 10 5 7 ; 
the nominal value implies = 10 7 . 
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TABLE 1 

Observational details for twenty-one new light curves of OGLE-TR-56b 



TV n ti t; l +■ 

_L 1 dllol L 


Tn Cltvi 1 TYlOTlt 
J_I±OLI LL11H_.11L 


± 1 Oj1111_.o 


1 JAN Ul) Lll 


Filter 


Ri nniTio' 

Lj 1111 1111 tL 


Read 


A 1 TTYl ace; 
nil iiicxoo 


Time observed 


Atm 


OL. Llllji, 


CUT) 






(sec) 






(sec) 




(before ingress) 


stability 


(arcsec) 


20030730r 


MaglC-SITe 


47 (0) 


- 


r' 


lxl 


23 


1.1-1.7 


2.4 hrs 


(+27 min) 


- 


- 


20030730z 


MaglC-SITe 


48 (0 


- 


z' 


1x1 


23 


1.1-1.7 


2.4 firs 


K -\-12 mm) 


- 


- 


20050419B 


MaglC-SITe 


59 (0) 


- 


B 


lxl 


23 


1.0-1.6 


o.z nrs 


, — 54 min) 


- 


- 


200504191 


MaglC-SITe 


59 (0 


— 




lxl 


23 


1.0-1.6 


3.3 hrs 


^ Q »yi 1 Ti \ 

i — OO 111111 t 


— 


- 


20050730 


MaglC-SITe 


140 (0) 


— 


1 


lxl 


23 


1.0-1.4 


3.2 hrs 


(-1-10 min) 


— 


— 


20050731 


MaglC-bf lo 


113 (0) 




1 


lxl 


23 


1.1-1.9 


3.0 hrs 


(+37 min) 






20060622 


POETS 


6642 (121) b 


2 


1s 


lxl 




1.0-1.4 


3.9 hrs 


-j-13 min) 


Stable 


0.25-0.5 


2006062b 


Maglu-blle 


245 (0) 






lxl 


23 


1.0—2.3 


4.6 hrs 


(+71 min) 






20060714 


POETS 


2200 (588) d 


1 


Is 


lxl 




1.0-1.4 


4.6 hrs 


(+88 min) 


Variable 


0.25-0.5 


200607201N 


FU-blb 


1867 (1379) 


6 


None 


lxl 




1.0—1.1 


5.5 hrs 


(+17 min) 


Stable 


0.4 


20070830 


IMACS 


352 (l) f ' 


15, 30 


WB§ 


lxl 


25 


1.0-1.5 


4.4 hrs 


(+112 min) 


Stable 


0.6 


20080514 


MaglC-e2v 


187 (0) 


60 




lxl 


5 


1.0-1.5 


3.6 hrs 


(+56 min) 


Superb 


0.3-0.4 


20080612 


MagIC-e2v 


604 (64) h 


17, 20 




lxl 


5 


1.0-1.2 


7.9 hrs 


(+264 min) 


Stable 


0.3-0.5 


20080727 


MagIC-e2v 


672 (0) 1 


10,12 




2x2 


5 


1.0-1.3 


4.3 hrs 


(+100 min) 


Stable 


0.7 


20090428 


CMOS 


309 (0) 


12-30 




2x2 


24 


1.0-1.5 


3.4 hrs 


(+32 min) 


Stable 


0.4-0.6 


20090504 


MagIC-e2v 


623 (2) 1 


23, 25 




2x2 


C 


1.0-1.2 


5.5 hrs 


(+90 min) 


Variable 


0.5-1.0 


20090510 


MagIC-e2v 


460 (2) f 


20, 30 




2x2 


C 


1.0-1.2 


3.7 hrs 


(+94 min) 


Stable 


0.7 


20090521 


MagIC-e2v 


626 (0) 


22-25 




2x2 


c 


1.0-1.2 


4.6 hrs 


(+65 min) 


Variable 


0.4-0.7 


20090612 


MagIC-e2v 


772 (3) f 


17, 20 




2x2 


c 


1.0-2.3 


4.3 hrs 


(+63 min) 


Stable 


0.4-0.7 


20090613 


MagIC-e2v 


715 (0) 


20 




2x2 


c 


1.0-2.0 


4.3 hrs 


(+48 min) 


Variable 


0.6-0.8 


20090728 


MagIC-e2v 


535 (5) f 


30-60 


i' 


2x2 


c 


1.0-1.2 


5.3 hrs 


(+187 min) 


Variable 


0.4-0.9 



a Number of frames used (additional frames that were discarded), 
k Discarded during meridian crossing (images elongated). 
c Readout is a few miliscconds in frame transfer mode. 



Discarded baseline after transit at a different pointing, which did not return to the same ratio level as before. 
° Discarded 31 frames during meridian crossing, 1347 frames at airmass > 1.1, and 1 frame with an aberrant ratio. 
^ Discarded a few anomalous frames, e.g. due to a seeing spike on one or more frames 
g WB6300-9500, a wide band filter from 630-950 nm. 

k Discarded initial frames with low counts on target (50 points) and a minor seeing spike (14 points). 

1 Discarded frames with large target diameters (poor seeing): >11.5 pixel (20080727) and >11.0 pixel (20090504) 



TABLE 2 

Aperture photometry parameters for each transit 



Transit 


Comp. 


Aperturc a 


Sky radius, width Slope removed 


Scatter 13 


CUT) 


Stars 


(pixels) 


(pixels) 


(mmag) 



20060622 


5 


10° 








1.0 


20060714 


2 


8.2 


30, 


10 




1.7 


20060720R 


6 


6.0 


50, 


5 




1.0 


20060720V 


5 


7.0 


10, 


5 




1.2 


20060720N 


2 


14 


70, 


5 


-3.6%, -9.6% d 


0.7 


20070830 


3 


7.6 


20, 


20 


-0.96% ° 


1.1 


20080514 


5 


14.5 C 








0.6 


20080612 


5 


12.2 


50, 


10 




0.9 


20080727 


3 


9.0 


25, 


30 


+0.16% f 


1.3 


20090428 


4 


4.0 


30, 


5 




1.0 


20090504 


2 


8.2 


20, 


30 


+1.48% e 


1.1 


20090510 


3 


7.8 


10, 


10 




1.0 


20090521 


5 


7.0 


30, 


10 




0.9 


20090612 


5 


8.2 


30, 


10 




0.8 


20090613 


3 


10.2 


25, 


10 




1.1 


20090728 


2 


8.2 


15, 


20 




1.2 



a Except as noted, the aperture is a circular radius around the star. 

k Scatter on out-of-transit flux, binned to 2 minutes. 

c Box half-width for alternate photometry method described in § |2.l| 

^ Trend removed for each star, individually, against airmass, in units of Z day 

e Trend removed against time, in units of flux day -1 . 

^ Trend removed against seeing in pixels, in units of flux pixels -1 . 
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TABLE 3 

Flux values for new transits of OGLE-TR-56B a 



Mid-exposure (JDjjtc) Mid-exposure (BJDtdb) Flux 



Error 



2453908.670845 
2453908.670868 
2453908.670891 
2453908.670914 
2453908.670937 



2453908.677416 
2453908.677439 
2453908.677462 
2453908.677485 
2453908.677508 



1.002203 0.001641 

1.007551 0.001641 

0.9947374 0.001641 

1.003771 0.001641 

0.9926688 0.001641 



a Full tabic available online. 



TABLE 4 

New system parameters for OGLE-TR-56 



Parameter 



Value 3 



Radius ratio, k 
Semimajor axis ratio, a/R t 
Inclination, i (deg) 
Eccentricity, e 

Linear limb-darkening coefficient, Uni 
Quadratic limb-darkening coefficient, 
Impact parameter, b 
Total transit duration, T14 (sec) 
Duration of full transit, T23 (sec) 
Ingress/egress duration, T12 = T34 (sec) 
Orbital period, P (days) 
Reference epoch, To (BJDtdb) 
Semimajor axis, a (AU) 
Planetary radius, R p (Rj) 
Planetary mass, M p (Mj) 
Stellar radius, R* (Rq) 
Stellar mass, M, (M Q ) 



0.1039 ±0.0004 
3.065 ±0.022 
73.72 ± 0.18 
(fixed) 
0.169 ±0.033 
0.175 ±0.041 
0.859 ± 0.003 
7931 ± 38 
2896±1° 0U 
251Btg 

1.21191096 ±0.00000065 
2453936.60070 ± 0.00044 
0.01942 ± 0.00015 b 
1.378 ±0.090 b 
1.39±0.18 b 
1.363 ± 0.089 b 
1.228 ±0.078 b 



a Median value of parameter distribution from joint fit to 23 light curves, 

with errors reported from the 68.3% credible interv al. 

b Stellar mass and radius values from Torres et al.|i2008|. 



TABLE 5 

Parameters for independent fits to each transit of OGLE-TR-56B 



Transit 


TAR a 


k 


t 


(deg 


) 




a/R, 




T14 






T23 


T12 = T34 


20030730r 


1 


7 


0.10414 ± 0.0072 


72 


43 ± 


2 


8 


2.836 ±0.32 


8688 


0± 


696 





3292.1 ± 1271.0 


2693.0 ± 895.0 


20030730z 


1 





0.1039 (fixed) 


74 


05 ± 


1 


2 


3.135 ±0.17 


7707 


3± 


298 





2737.2 ± 637.0 


2489.3 ± 372.0 


20050730r+z b 


1 


6 


0.09856 ± 0.004 


74 


47 ± 


1 


9 


3.199 ±0.28 


7500 


7± 


417 





3044.2 ± 836.0 


2229.3 ± 558.0 


20050731 


1 


3 


0.10653 ±0.0091 


74 


21 ± 


1 


7 


3.156 ±0.21 


7728 


1 ± 


330 





2749.0 ± 804.0 


2485.9 ± 499.0 


20060622 


1 


9 


0.10392 ±0.0023 


74 


15 ± 


1 


1 


3.077 ± 0.14 


8143 


4± 


216 





3510.9 ± 516.0 


2317.3 ±334.0 


20060626 


2 


1 


0.09986 ± 0.0064 


74 


12 ± 


3 


9 


3.095 ± 0.54 


7938 


0± 


696 





3446.9 ± 1528.0 


2250.2 ± 1046.0 


20060720P C 


1 





0.10269 ±0.0017 


77 


48 ± 





8 


3.754 ± 0.15 


6817 


9± 


129 





3452.3 ± 281.0 


1684.0 ± 192.0 


20060720R 


1 


1 


0.10064 ± 0.0097 


73 


85 ± 


1 


9 


3.094 ±0.24 


7764 


9± 


361 





3001.0 ± 783.0 


2362.6 ±492.0 


20060720V 


1 


1 


0.1011 ± 0.0066 


76 


74 ± 


4 





3.503 ± 0.7 


7428 


3± 


475 





3973.5 ± 948.0 


1733.8 ± 664.0 


20060720V±R d 


1 


3 


0.1039 (fixed) 


73 


82 ± 





9 


3.092 ±0.14 


7822 


6± 


265 





2766.6 ± 449.0 


2528.1 ± 283.0 


20060720N 


1 


6 


0.1033 ±0.0021 


71 


12 ± 


1 


1 


3.116 ±0.15 


7858 


5± 


204 





3088.4 ± 604.0 


2385.5 ± 380.0 


20070830 


1 





0.10228 ±0.0018 


73 


91 ± 





9 


3.065 ±0.11 


8025 


8± 


182 





3290.0 ± 439.0 


2369.0 ± 283.0 


20080514 


1 


1 


0.10408 ±0.0012 


74 


17 ± 





5 


3.131 ± 0.07 


7817 


2± 


121 





3008.5 ± 276.0 


2405.0 ± 180.0 


20080612 


1 


8 


0.10419 ±0.0029 


73 


66 ± 


1 


1 


3.061 ± 0.18 


7924 


2± 


293 





2807.6 ± 857.0 


2556.6 ±517.0 


20080727 


1 


7 


0.10695 ±0.0045 


73 


84 ± 


1 


6 


3.08 ±0.2 


7968 


5± 


346 





2831.0 ± 961.0 


2563.7 ± 576.0 


20090428 


1 


7 


0.10751 ± 0.0108 


72 


74 ± 


1 


6 


2.987 ± 0.17 


7829 


0± 


296 





1676.7 ± 935.0 


3051.5 ±528.0 


20090504 


2 


6 


0.10244 ± 0.0042 


75 


23 ± 


2 


3 


3.243 ± 0.33 


7817 


5± 


404 





3712.0 ± 905.0 


2053.5 ±605.0 


20090510 


2 





0.10861 ± 0.0037 


73 


46 ± 


1 


5 


3.013 ±0.19 


8182 


1 ± 


340 





2795.2 ± 1012.0 


2691.6 ± 630.0 


20090521 


1 


3 


0.10159 ±0.0018 


74 


01 ± 





8 


3.085 ±0.1 


7954 


7± 


171 





3286.5 ±391.0 


2335.0 ± 253.0 


20090612 


1 


9 


0.10597 ± 0.0025 


73 


88 ± 


1 


1 


3.07 ± 0.15 


8039 


0± 


245 





3048.9 ± 598.0 


2496.1 ± 387.0 


20090613 


1 


9 


0.10429 ±0.0039 


71 


67 ± 


2 


2 


3.188 ±0.32 


7801 


5± 


433 





3291.9 ± 926.0 


2258.6 ±641.0 


20090728 


2 


1 


0.09963 ± 0.0039 


75 


09 ± 


2 


2 


3.271 ± 0.33 


7524 


5± 


408 





3377.5 ± 857.0 


2076.4 ± 585.0 



a Multiplicative factor indicating excess noise, as derived from the time-averaged residual method, 
k Combined r' and z' photometry. 

c New fit to original photometry of |Pont et al.| | |2b07) 
Combined V and R photometry. 
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TABLE 6 

OGLE-TR-56B Transit Midtimes and Residuals 



Transit 


Number 


Midtime (BJD TDB ) 


yj-Kj (s ) 


a 


Source 


20010615 a 


-1536 


2452075 


.10534 ±0.00170 


-11 ± 147 


-0.1 


|Torres et al.|(|2004|| 


20030730r 


-896 


2452850 


.73111 ± 0.00092 


226 ± 80 


2.9 


this work 


20030730z 


-896 


2452850 


.72730 ± 0.00083 


-103 ± 72 


-1.4 


this work 


20030730r+z b 


-896 


2452850 


.72908 ± 0.00045 


51 ± 39 


1. 1 


this work 


20050419B 


-377 


2453479 


.71055 ± 0.00540 


24 ± 466 


0.1 


this work 


200504191 


-377 


2453479 


.70867 ± 0.00290 


— 139 ± 246 


-0.6 


this work 


20050730 


-293 


2453581 


.51076 ± 0.00096 


-3 ± 83 


-0.0 


this work 


20050731 


-292 


2453582 


.72328 ± 0.00073 


49 ± 63 


0.8 


this work 


20060622 


-23 


2453908 


.72603 ± 0.00060 


-63 ± 51 


-1.2 


this work 


20060626 


-19 


2453913 


.57171 ± 0.00130 


-233 ± 114 


-2.0 


this work 


20060714 


-5 


2453930 


.54256 ± 0.00170 


121 ± 148 


0.8 


this work 


20060720P C 





2453936 


.59875 ± 0.00100 


-169 ± 86 


-1.9 


|Pont et al.| <|2007J 


20060720V 





2453936 


.59735 ± 0.00080 


-291 ± 69 


-4.2 


this work 


20060720R 





2453936 


.60226 ± 0.00072 


134 ± 62 


2.2 


this work 


zuuuu i £\i v T^rx 


o 


2453936 


.59973 ± 0.00092 


—85 ± 79 


-1.1 


Llllo WU1 rV 




o 


2453936 


.59908 ± 0.00045 


— 141 ± 39 


-3.6 


tnio ~\iiTf~\r\r 
bLlla WUI rv 


20070830 


335 


2454342 


.59179 ± 0.00046 


78 ±40 


2.0 


this work 


20080514 


548 


2454600 


.72852 ±0.00031 


52 ±27 


1.9 


this work 


20080612 


572 


2454629 


.81483 ± 0.00083 


90 ± 71 


1.3 


this work 


20080727 


609 


2454674 


.65443 ± 0.00085 


-5± 73 


-0.1 


this work 


20090428 


836 


2454949 


.75803 ± 0.00065 


-21 ± 56 


-0.4 


this work 


20090504 


841 


2454955 


.81712 ±0.00110 


-61 ± 92 


-0.7 


this work 


20090510 


846 


2454961 


.87779 ± 0.00069 


35 ±60 


0.6 


this work 


20090521 


855 


2454972 


.78553 ± 0.00044 


81 ± 38 


2.2 


this work 


20090612 


873 


2454994 


.59691 ± 0.00060 


-179 ±52 


-3.4 


this work 


20090613 


874 


2454995 


.80961 ± 0.00079 


-111 ±68 


-1.6 


this work 


20090728 


911 


2455040 


.65158 ±0.00093 


-1 ± 80 


-0.0 


this work 



a Time has been adjusted from the published times into the BJDtt time system by adding UTC-TT 
conversion of 64.184 sec. 

Midtime value w hen the r' and z' light curves are combined and fit independently. 
c Light curve from |Pont et al~| ) |2007| ; the midtime reported is from our independent fit of the original 
photometry. 

Midtime value when the R and V light curves are combined and fit independently. 
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Fig. 1. — Twenty-three light curves observed over nineteen different transit epochs of OGLE-TR-56b. All data are binned to 2 minutes 
to aid comparison. The solid lines show the best model for all curves, color-coded by the instrument and filter used: SITe/B = black, 
FORSl/V=purple, SITe/r' = dark blue, FORSl/K=light blue, IMACS/WB6500-9300 = cyan, e2v/i' = green, GMOS/i' = dark green, 
POETS/is = yellow, SITe/7 = orange, SITe/z' = red, and POETS/no filter = gray. The residuals from the model are plotted below each 
curve. 
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Fig. 2. — The field of view for OGLE-TR-56, as observed on 2008 May 14 with MagIC-e2v on the 6.5 m Baade telescope. Field of view 
is 38" x 38", with north up and east to the left. 
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Fig. 3. — Parameter variation of individual transits of OGLE-TR-56b. Data from the individual MCMC fits (Table [5jl is shown 
for the inclination (i), semimajor axis (a/R t ), planet-star radius ratio (fc), total transit duration (T14), full transit duration (T23), and 
ingress/egress duration, assuming a circular orbit (T12 = T34). Note that all errors have been scaled upward based on the factor calculated 
from residual permutation. For comparison, the value derived from th e joint MCMC fit is plotted as a solid blue lines with dashed ±1<t 
errors. We highlight our independent fit to the light curve published by |Pont et al,| ( |2007} in red to illustrate that it is markedly discrepant 
for all parameters except the radius ratio. 
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Fig. 4. — Three light curves of the transit on 2006 July 20, observed independently by two instruments, with a fourth curve on a different 
night for comparison. Data has not been binned. Top curve: original image-subtraction photometry for transit 20060720P observed with 
FORS1 on the VLT by |Pont et al.| ( |2007| l, with _R-band data in red and V-band data in black. Second curve: the same VLT transit, redone 
with aperture photometry fo r thi s work. Third curve: independent observation of the same transit with no filter, 20060720N, using POETS 
on Magellan (discussed in § |2.1| . The bottom curve is a different transit, 20080514, to illustrate the agreement of all but the top curve 
(20067020P) w ith the othe r transits presented in this work. Two models are shown in each plot: a gray dashed line for the best fit to the 
photometry by |Pont et al,||2007[ l (refit for this work), and a black solid line showing the best joint-fit to all 23 new light curves. 
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Orbital Periods 

Fig. 5. — Observed minus calculated midtimes for OGLE-TR-56b. Timing residuals using the new ephemeris (Equation^. Solid symbols 
mark transits that were used to calcul ate the ephemeris . The solid red triangle at —1536 is the OGLE midtime, while the open red triangle 
at is our fit to the photometry from |Pont et al.| ( |2007} . The two solid blue squares at -896 and represent independent fits to combined- 
filter light curves, which were used in place ot the pair of open circles at those same values, representing each individual filter curve. See 
§|4]for more details. The bottom panel shows a zoomed-in view of the region around 0. The solid line shows the expected time of transit, 
wmle the dashed lines represent the la and 3a errors on the calculated orbital period, indicating the slopes that result for a mis-determined 
period. Notice that the la and 3a error lines converge around the zero epoch and therefore appear practically indistinguishable in the 
lower plot. 



